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ABSTRACT 

We measure the mass and size of cloud fragments in several molecular clouds continuously over a wide 
range of spatial scales (0.05 < r/pc < 3). Based on the recently developed "dendrogram-technique", 
this characterizes dense cores as well as the enveloping clouds. "Larson's 3''^ Law" of constant column 
density, m(r) oc r^, is not well suited to describe the derived mass-size data. Solar neighborhood 
clouds not forming massive stars (< 10 M©; Pipe Nebula, Taurus, Perseus, and Ophiuchus) obey 



m(r) < 870 Af0 (r/pc) 



1.33 



In contrast to this, clouds forming massive stars (Orion A, GIO. 15—0.34, Gil. 11— 0.12) do exceed the 
aforementioned relation. Thus, this limiting mass-size relation may approximate a threshold for the 
formation of massive stars. Across all clouds, cluster-forming cloud fragments are found to be — at 
given radius — more massive than fragments devoid of clusters. The cluster-bearing fragments are 
found to roughly obey a mass-size law m oc r^-'^'^ (where the exponent is highly uncertain in any given 
cloud, but is certainly smaller than 1.5). 

Subject headings: ISM: clouds; methods: data analysis; stars: formation 



1. INTRODUCTION 

Most of our present understanding of star formation 
processes is based on detailed studies of solar neighbor- 
hood molecular clouds (closer ~ 500 pc). To this end 
past research has, e.g., studied the masses and sizes of 
dense cores in molecular clouds (< 0.1 pc size) such 
as Perseus, Taurus, Oph i uchus , Orion, and the Pipe 
Nebula (e.g. |M otte et al.' '1998', Johnst one et all 120001 
iHatcheU et "all 2005. Enoch ct al. 2007). Further re- 
search studied clumps in these clouds (some 0.1 pc) 
and the clouds (> 10 pc) containing the cores (e.g. , 
IWilliams etUlfTaM ICambresvl [19991 iKirk et all [20061 : 



Most subsequent work refers to this relation as "Larson's 
3'''^ la w" , and replaces the orig inal result with m(r) oc 
(e.g., iMcKee fc Ostrikeii [20071 ) . This "law of constant 
column density" (with respect to scale, r) is now con- 
sidered one of the fundamental propert ies of molecu- 
lar cloud struct ure (e.g.. Ballesteros-Paredes et al.ll2007l 
IMcKee fc Os~ik cr 2007. iBergin fc Tafallall2007ir We do, 
however, not know whether this relation is still consis- 
tent with up-to-date column density maps of molecular 
clouds. 



see 

and clouds). This research does, however, not probe the 
relation between the properties of cores, clumps, and 
clouds: traditionally, every domain is characterized and 
analyzed separately. As a result, even in the solar neigh- 
borhood, it is still not known how the core densities (and 
thus star-formation properties) relate to the state of the 
surrounding cloud. 

To be precise, we do in principle know a bit about 
the relation between cloud structure at large an d small 
scale. For structure within molecular clouds, iLarsonl 
(jl981|) concluded (in his Eq. 5) that the mass contained 
within the radius r obeys a power-law. 



Part I of the present series (jKauffmann et all I2010D 
describes a new technique to extract mass-size relations 



Wilhams et al.||2000| for definitions of cores, clumps, from cloud maps. It is based on "dendrograms 



m(r) = 460 Mq (r/pc) 



1.9 



(1) 
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, a tree- 
based segmentation of cloud structure ( R osolowskv et al] 
|2008[) . Here, we employ this technique to study the 
molecular clouds in Perseus, Taurus, Ophiuchus, Orion, 
and the Pipe Nebula. To illustrate the properties of more 
massive clouds, we also include data for two more distant 
clouds of high density (farther than 2 kpc; GIO. 15—0.34 
and Gil. 11-0. 12). 

The present paper summarizes the analysis method in 
Sec. [21 The main quantitative analysis of the maps is pre- 
sented in Sec. [3] Section |4| systematizes the results and 
interprets them in the context of our present knowledge 
of star formation regions. This discussion is supported by 
model calculations in Appendices \K\ and [Bl We conclude 
with a summary in Sec. [5] 

2. METHOD & DATA 

2.1. Data Processing 

2.1.1. Basic Map Analysis 

Our basic analysis approach is summarized in Sec. 2.1 
of part I, and illustrated in Fig. 1 of the same paper. In 
essence, starting from a set of local maxima, we contour 
a given column density map at all levels possible. For 
every contour, we derive the enclosed mass and area, A. 
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Fig. 1. — Mass-size relations for the clouds listed in Table [T] In most clouds, data from two different observational techniques (e.g., 
dust extinction and emission) are combined to provide a comprehensive picture probing a wide range of spatial scales. Solid black lines 
highlight cloud regions containing the most massiv e cloud fragment found for small radii. Green lines starting in circles ind icate the 
global mass-size relations discussed in Section 12.1.21 The other dotted lines give reference mass-size relations as discussed in Sec. I2.1."3l In 
Orion A, published data extracted using CLUMPFIND-like approaches are plotted [diamonds mark SCUBA data, crosses indicate CO 
observations) instead of using the contour-based scheme employed here. Further, the triangle indicates the large-scale Orion A extinction 
mass measurement discussed in the text. 



Following the terminology of lPeretto fc FulleJ (|2009[ ). we 
define "cloud fragments" in the maps as such regions 
enclosed by a continuous column density contour. The 
area is used to derive effective radii, 

r = (A/7r)i/2 . (2) 

Subsequent contours are usually nested in the map. This 
defines a relation between measurements. This essen- 
tially yields series of mass-size measurements. In the 
plots shown in this paper, such series are drawn using 
continuous lines (e.g., Fig.[T]). 

In practice, the processing is implemented using 
the dendr ogram m ap analysis technique introduced by 
[Ro solows kv et al.l ([2008). As a bonus, this also yields 
diagnostic diagrams on the cloud hierarchy. We do not 
use this feature in the present paper, though. 

In our analysis, we reject all cloud fragments that have 
a diameter (i.e., 2r) smaller than the map resolution. 
Further, we require a minimum contrast between the lo- 



cal maxima used to seed the contouring. Here, this limit 
is set to the noise level times a factor 3. Further, we do 
not characterize column densities below a certain thresh- 
old. As discussed below. Table [1] lists these parameters 
for each map studied here. 

2.1.2. Calculation of Mass-Size Slopes 

Part I demonstrates that mass-size data can conve- 
niently be characterized using power-laws of the form 

m{r) = mo (r/pc)'' . (3) 

These are characterized by slope, 6, and intercept, niQ. 
As shown in part I, slopes can be derived using various 
methods. As demonstrated in part I, 5 < 2 in our work, 
since b > 2 would imply an increase of the mean column 
density with radius. 

We define global slopes to capture trends between 
small and large spatial scales (Sec. 4.3 in part I). To cal- 
culate these, we derive the maximum fragment masses. 
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"T-max, observed at radii of rgm = 0.05 pc and rig = 5.0 pc. 
(These radii are chosen to permit comparison between 
clouds, as becomes more obvious below.) Based on these, 
we derive the global slope. 



A data summary for our sample is provided in Table [TJ 
This includes parameters used for the source extraction. 



bo-ioh — 



In [TOmax (rig ) /mmax (r 

ln[rig/rsin] 



(4) 



As illustrated in Fig. [U this slope is defined such that 
m(r) cx r^siob connects the mass and size measurements 
at 0.05 pc and 5.0 pc radius (for appropriate intercept). 

To characterize trends at a given spatial scale, we use 
tangential slopes (Sec. 4.4 in part I). These are derived 
infinitesimally at a given radius, 



b{r) 



d In(m[r']) 



d In(r') 



(5) 



As illustrated in Fig. 7 of part I, the measured slopes are 
smoothed and filtered to improve the data quality. 

2.1.3. Reference Mass-Size Relations 

Section 3.1 of part I introduces various reference mass- 
size relations. These can be used to navigate more intu- 
itively within the data space. They are thus featured in 
most figures of the present paper (e.g.. Fig. [T]). Specifi- 
cally, lines of constant mean column density obey a mass- 
size law m{r) oc . Our plots contain relations for mean 
column densities separated by factors of 10. Static equi- 
librium models of isothermal spheres, on the other hand, 
obey m(r) oc r. We draw a model relation for a gas 
temperature of 10 K. 

We stress that we obtain two-dimensional mass-size 
relations from column density maps. These are related 
to, but not identical with, mass-size laws obtained from 
three-dimensional density maps. This is further explored 
in Sec. EJl 



2.2. Combining Data for several Clouds 

Part I describes in detail how dust extinction and emis- 
sion maps can be used to derive column density maps for 
molecular cloud complexes. Once calibrated to a com- 
mon mass conversion scale (Sec. 4.2 of part I), dust ex- 
tinction and emission data for a given cloud can be com- 
bined to probe the masses and sizes across a vast range 
of spatial scales. In part I, we did this for Perseus only. 
Here, we combine our Perseus data with maps of other 
clouds. We do this by repeating the analysis already 
carried our for Perseus. 

We start with a survey of a sample of well-studied 
nearby clouds that are essentially devoid of high-mass 
stars. Besides the Perseus molecular cloud, those in 
Taurus, Ophiuchus, and the Pipe Nebula are examined 
here. As a general reference to clouds also forming high- 
mass stars, we include the Orion A cloud. Then, the 
more remote GIO. 15-0.34 (hereafter GIO; - 2.1 kpc) 
and Gil. 11— 0.12 complexes (hereafter Gil; ^ 3.6 kpc) 
are studied to build a first bridge towards the study 
of r elatively distant s ites of high mass star formation 
(see lPillai et al.l (|2007f) for distances and references; GIO 
is further discussed by Wood & Churchwell (1989) and 
iThompson et al.l 120061 while .Pillai et al.. (,2006b. ) study 
Gi l). Note that Gil is an I nfrared Dark Cloud (IRDC; 
see iMenten et all 120051 and iBeuther etldl 120071 for re- 
views). 



2.2.1. Data and Analysis 

The extinction map analysis for Ophiuchus 
(iRidge et al.ii2006il and the Pipe Nebula (jLombardi et al.l 
120061 ) is analogous to the one for the Perseus cloud 
carried out in part I. The visual extinction in the 
Pipe Nebula map is reduced by 1.34 mag (reduc- 
tion of 0.15 mag in Ak), following the analysis by 
fLombardi et al. ( 20081) . For T aurus, we use the map 
by iRowles fc FroebrichI (|2009l ). Its resolution varies 
throughout the map. This is analogous to a region- 
dependent smoothing, which may affect mass estimates 
(Sec. 4.1 in p art I). We thus use this map with caution. 
In Taurus ( Kauffmann et al.l |200^ and Ophiuchus 
(jEnoch et al.ll2007[ ) the processing of the dust emission 
maps follows the Perseus Bolocam analysis of part I. 
Note that the Taurus MAMBO maps do not cover all 
of the cloud and may therefore give a biased view of the 
dense core conditions. Pipe Nebula regions with high 
column density must be probed in detail. To do this, we 
include a SCUBA map for B68 (Alvcsct al. 2 00l, which 
is kiii dly provided by J. Alves. [Roman- Ztiniga et al.l 
(|2009f ) mapped the Pipe's B59 region in extinction. We 
process these maps using our analysis scheme. For all 
sources. Table [1] lists the free parameters used in our 
source extraction algorithm. 

The resulting mass-size relations form part of Fig. [TJ 
To illustrate some aspects of the spatial cloud structure. 
Fig. [H highlights (in black) for every cloud the mass-size 
evolution of the most massive fragment found at small 
radii. Note that, in Taurus and the Pipe Nebula, these 
fragments do not form part of the largest cloud fragments 
found in the maps. 

In Orion A, no data suited for a reliable dendrogram 
structure analysis o f column densities is yet a vailable. 
Instead, we use the IRowles k FroebrichI (|2009l ) extinc- 
tion maps to derive a single mass and size measure- 
ment on the largest scales probed by that map. On 
smaller scales, the extinction map is not reliable because 
of too few background sources. For reference, we also 
plot mass and size measureme nts for ^'^CO clum ps in 
Orion, as published elsewhere (jBallv et al.l Il987t their 
Table 1, plus text statements). These data are, how- 
ever, not used in the quantitative analysis. At smaller 
scales, we fold in published data from SCUBA stud- 
ies employing CLUMPFIND-like data analysis methods 
(|Nutter fc Ward-ThomT)Sonll2007[ ). As illustrated in Fig. 
4 of part I, at given radius, these set a lower limit to the 
maximum dendrogram-derived masses usually employed 
here. 

For GIO (the data are kindly provided — in advance 
of publication — by M. Thompson, J. Hatchell, and F. 
Wyrowski) and Gil (ICarev et al.ll2000D . we run our con- 
tour analysis on SCUBA maps to derive mass-size data 
on large spatial scales. In Gil, the structure on very 
small scales is probed in a similar fashion, using inter- 
ferometric dust emission maps obtained using the Sub- 
millimeter Array (SMA; Pillai et al., in prep.). Unfor- 
tunately, the SMA map does not cover all of Gil, and 
the derived mass-size data are likely to be biased (i.e., 
the most dense core is not covered by our data). We do 
not use our existing SMA data for GIO, since the data is 
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TABLE 1 

Data included in this study. 



Region / Distance 


Data 


Resolution 


Noise 
mag 


Rejection Threshold 
mag 


Reference 


clouds only forming low-mass stars: 










Taurus 


extinction 


< 6' 


0.4 


2.0 


Rowles & Froebrich (■20091 


d = 140 pc 


dust emission (10.0 K) 


20" 


1.4 


4.3 


Kauffmann et al. C2008') 


Perseus 


extinction 


5' 


0.4 


2.0 


Ridge et al. (2006) 


d = 260 pc 


dust emission (12.5 K) 


30" 


1.4 


4.3 


Enoch et al. (20061 


Ophiuchus 


extinction 





0.6 


2.0 


Kiage et ai. (zUUb) 


d = 120 pc 


dust emission (12.5 K) 


30" 


2.5 


7.6 




Pipe Nebula 


extinction 


dU 


U.o 


4.U 


LjomDarai et ai. izuuo) 


d = 130 pc 












B59 


extinction 


zu 


1 2 


o.u 


XXUIlidll-ZjLlIiliiti cli cLi. I ZUUy ) 


d = 130 pc 










B68 


dust emission (10.0 K) 


15" 


1.1 


3.3 


J. Alves, priv. comm. 


d = 130 pc 












clouds also forming high-mass stars: 










Orion A 


"CO 


1.'7 






Ballv et al. (1987) 


d = 414 pc 


dust emission (20.0 K) 


14" 






Nutter & Ward-Thompson f2007) 


GlO.15-0.34 


dust emission (20.0 K) 


14" 


6.4 


19.2 


Thompson et ah, in prep. 


d = 2100 pc 












Gil. 11-0.12 


dust emission (15.0 K) 


14" 


6.6 


19.7 


Carev et ah f2000^ 


d = 3600 pc 


dust emission (15.0 K) 


4" 


4.9 


14.6 


Pillai et ah, in prep. 



misleading in this context (our GIO map only covers a mi- 
nor dust emission peak). In the SCUBA and SMA maps, 
structure larger than a certain spatial scale is removed 
during data reduction. Setting other uncertainties aside, 
the true intensities and column densities will therefore 
be larger than derived here. The adopted dust te mpera- 
tures are inspired bv iPillai et al.i (,2006a) and .PiUai et alJ 
(|200l . 

2.2.2. Mass Estimates 

Table [T] lists (as part of the data column) the dust 
temperatures adopted for each cloud. The temperatures 
for individual fragments may deviate from these mean 
values by several Kelvin. This introduces an associated 
uncertainty in mass estimates from dust emission. An 
uncertainty of order 25% is probably reasonable for ex- 
treme cases. 

For Taurus, Figure [T] reveals a jump of order of a factor 
2 in the mass-size relation from dust extinction to dust 
emission data, with the dust emission data implying the 
higher masses. This suggests a problem in the relative 
calibration of tracers, or a problem with the generation of 
the extinction ma]:3 Here, we use the Taurus extinction 
map with caution in our analysis. 

B59 of the Pipe Nebula map appears to suffer from 
even larger mass biase s; for ex actly the same region, 
when compared to the IRoman-Z iifiiga ct al. (200!|) re- 
sults, the lower-resolution iLombardi et al.l (|2006i) map 

^ The other extinction maps used here — created by Lada, Alves, 
Lombardi, and collaborators — use a spatially constant resolution. 
Also, the calculation of the extinction for a given star differs in 
details. 



implies masses lower by a factor ~ 3. Given the extreme 
column densities in the target region (B59), this result is 
not entirely surprising. We do not expect such biases in 
the other regions (where areas of high column density are 
anyway probed by dust emission). Still, the B59 results 
suggest to use extinction maps with extreme caution. 

3. MASS-SIZE RELATIONS: A CLOUD SAMPLE 

At this point, all mass-size data has been collected and 
processed. In the following sections, we highlight three 
particular trends seen in this sample data. 

In principle, the mass-size data for a given cloud might 
significantly depend on the viewing direction. Since we 
can only observe the cloud projections as seen from earth, 
our data might be biased because of our specific viewing 
direction. Here, we make the assumption that this is 
not the case. This notion is supported by our obser- 
vational findings: all nearby clouds not forming mas- 
sive stars have (in a broad sense) similar properties. 
Also, mass-size differences between clouds correlate with 
projection-independent cloud properties (e.g., for given 
radius, clouds forming massive stars are more massive). 
This suggests that projection effects do not significantly 
bias the mass-size laws derived here. 

3.1. A limiting Mass-Size Relation for Massive Star 
Formation? 

As a first characterization of the cloud sample, we look 
at the maximum radius-dependent mass of cloud frag- 
ments. We begin with clouds not forming massive stars. 
In the 0.01 < r/pc < 10 radius range, essentially all 
fragments in such clouds have a mass smaller than some 
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limiting law, 

m(r) = 870 Mq (r/pc)!-^^ . (6) 

Figure [5] gives an illustration. This limit thus gives the 
typical mass range for structure in clouds like Taurus, 
Perseus, Ophiuchus, and the Pipe Nebula. 

In detail, this relation excludes a bright fragment in 
NGC1333, which we remove from our analysis as a pos- 
sibly unphysical outlier (its estimated mass exceeds those 
of other fragments of similar size by a factor 2, possibly 
because of neglected protostellar heating). The limit- 
ing law is derived by searching for the smallest intercept 
for which niQ (r/pc)** does still provide an upper limit 
to the data. Practically, this is done by varying b until 
max[m(r)/(r/pc)^] is minimized. 

Interestingly, our sample clouds with massive star for- 
mation exceed the limiting relation (Eq. iB]). This is 
shown in Fig. Wlb). For the Orion A cloud and GIO, 
we derive an excess of up to a factor 10 in the 0.01 < 
r/pc < 2 radius range. This suggests that these clouds 
have a structure significantly different from what is found 
for clouds not containing such stars (Fig. [Ha]). In this 
light, Eq. ^ may approximate a limit for massive star 
formation: it could be that only clouds also containing 
fragments that exceed Eq. ((6]) are able to form massive 
stars. Larger samples of clouds forming massive stars 
must be screened to prove this point. Note, however, that 
most fragments in Orion A do fulfill Eq. ©. This cloud 
does therefore also contain objects that have masses and 
sizes not distinguishable from those found for clouds not 
forming massive stars. 

As a mass-size limit for clouds not forming massive 
stars, Eq. (jH) is probably uncertain to just a few 10%. If 
we, e.g., not use Ophiuchus data in the derivation of Eq. 
([6]), then Ophiuchus would exceed the resulting mass-size 
limit by 30%. If we do the same with Perseus data, the 
excess is 15%. It is plausible to expect similar changes 
for other regions not forming massive stars. Observa- 
tional uncertainties are most likely of a similar order, if 
one adopts the mass measurement techniques used here 
(i.e., dust emission and extinction). Then, several un- 
certainties (e.g., dust opacities) are simply removed by 
calibration to the same standard. In an absolute sense 
(i.e., when considering the true masses), Eq. ([6]) is as ac- 
curate as the mass conversion standards used here (e.g., 
relation of dust emission and mass) . These are probably 
uncertain to less than a factor 2. 

We stress that, excluding observational uncertainties, 
Eq. Q provides a strict upper mass limit to the solar 
neighborhood sample as we have defined it here. Clouds 
violating Eq. ([6]) are not similar to the clouds in the solar 
neighborhood sample discussed here. This statement is 
sufficient for many purposes. 

3.2. Local Clouds obey similar Mass-Size Relations at 
large Scales 

Inspection of Figs.[T]and[5]suggests that all solar neigh- 
borhood clouds have similar masses at given size for radii 
> 1 pc. Specifically, the most massive fragments at given 
radius are essentially all within ±30% of 

m(r) = 400 Mq (r/pc)^-'^ (7) 
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Effective Radius [pc] 

Fig. 2. — Joint mass-size plot for all clouds, separated into clouds 
with (bottom) and without (top) active formation of high mass 
stars. See Fig. [l] for color coding. The shading indicates the ap- 
proximate mass limit for fragments in clouds without high mass 
star formation (top panel; Eq.|6]l. In our sample, all clouds forming 
high ma ss st ars contain fragments that exceed this limit (bottom). 
See Sec. l3.ll for details. 

in the 1 < r/pc < 4 radius range, when considering 
Taurus, Ophiuchus, Perseus, and the Pipe Nebula. Also 
our preliminary Orion A data is, at 11 pc radius, within 
40% of this law — with the caveat that we cannot exam- 
ine whether Orion A follows Eq. ([7]) down to 1 pc 
radius. It thus appears that — with Orion A as a possible 
exception — all local clouds are similar in their large-scale 
mass structure. 

This mass-siz e relation is ve ry similar to the one origi- 
nally derived bv lLarsonI (|1981f ). m(r) = 460 (r/pc)^-^ 
(see Eq. [I] above) . Our study thus confirms his result — 
but only for spatial scales ^ 1 pc. As we show through- 
out this paper (e.g., Eq. [6] and Figs. |4] and [5]), no single 
mass-size relation describes all structural aspects of our 
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observational data. 

Some cloud fragments ^ 1 pc are, however, much more 
massive than solar neighborhood clouds of similar size. 
GIO, for example, violates Eq. ([7]) by about an order of 
magnitude. The similarity between local clouds implied 
by Eq. (O might thus only hold for the solar neighbor- 
hood. 

3.3. Clusters vs. Isolated Star Formation 

It is obvious that the formation of a cluster requires 
a larger mass reservoir than necessary to form a single 
isolated star. Thus, one might naively expect that re- 
gions containing clusters are more massive than those 
devoid of stellar groups. If true, one would thus expect 
that, within a given cluster- forming cloud, the regions 
containing clusters are more massive than cluster-less re- 
gions of similar size. This hypothesis is tested in Sec. 
13.3.11 Also, one would expect that cluster-forming cloud 
fragments are more massive than all similar-sized frag- 
ments in clouds not containing clusters. Section 13.3.21 
investigates this issue. 

In our sample, the Pipe Nebula and Taurus serve as 
examples of regions dominated by isolated star forma- 
tion. Actually, except f or B59, the Pipe Ne bula does 
hardly form stars at all ([Forbrich et al.ll2009D . Perseus 
and Ophiuchus serve as examples for cluster-forming re- 
gions. They do contain clusters much more significant 
than any stellar aggregate found in Taurus and the Pipe 
NebulcO- Orion A is another example of a cluster-forming 
cloud. 

3.3.1. Cluster-forming Fragments dominate their Host 

Cloud 

We examine whether cluster-forming cloud fragments 
dominate the mass reservoir of their host cloud at all 
radii. This is executed in Fig. [TJ Here, we only con- 
sider cluster-forming clouds with high quality data, i.e. 
Perseus and Ophiuchus. Inspection of the column density 
maps reveals that the (highlighted) most massive small- 
scale features in Perseus and Ophiuchus are located in 
the NGC1333 and L1688 clusters, respectively. At small 
spatial scales, the most massive fragments are thus in- 
deed located in cluster-forming regions. Note, though, 
that cluster-forming regions do also contain fragments of 
lower mass. This is illustrated in the rightmost panel of 
Fig. 6 in part I: Bolocam-detected fragments in NGC1333 
cover a wide range in mass. 

In Ophiuchus, when examining larger spatial scales, 
fragments containing L1688 continue to be the most mas- 
sive ones. This is not exactly true for NGC1333 in 
Perseus, though. In this cloud, when considering ex- 
tinction maps, other fragments (drawn in red) are — by 
a small margin — the most massive ones at given radius. 
However, closer inspection (not presented here in detail) 

^ Clusters in near by clouds were recently surveyed by 
IGutermuth et all l|2009l ) . Based on Spitzer data, Ophiuchus and 
Perseus are found to contain clusters with > 130 members. The 
same study gives ~ 40 members for L1495 in Taurus. This latter 
value is in line with previous studies and s ets an upper limit to the 
size of stellar groups in Taurus (Table 4 of lKenvon eralll2008l ') . In 
the Pipe Nebula, the star formation activity is dominated by the 
B59 region (jF orbrich ct al. 2009). This group has ~ 20 members 
IjErooke et al.ll2007i ) . Taurus and the Pipe Nebula do thus not con- 
tain clusters as significant as the ones in Perseus and Ophiuchus. 
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Fig. 3. — Perseus and Ophiuchus cloud regions containing the 
most massive cloud fragment found for small radii (highlighted by 
dotted black lines) in comparison to structure in T auru s and the 
Pipe Nebula. See Fig. [T] for color coding. Section 13.31 describes 
that cluster-forming fragments exceed, at given size, those without 
clusters in mass. 

reveals that these fragments all contain the cluster IC348. 
In summary, cluster-forming fragments do thus indeed 
constitute the most massive cloud features at given ra- 
dius. Some cluster-bearing fragment, which constitutes 
the most massive cloud fragment at some radius, might 
however at a different radius be less massive than some 
other cluster-forming fragment. In this context, note 
that IC348 is probably much older than NGC1333 (e.g., 
'Gute rmuth et al1l2009t ). Dense star- forming cores, which 
manifest as small objects of large mass, are thus actually 
not expected to remain in the IC348 region. 

3.3.2. Cluster-forming Fragments exceed Fragments 
without Clusters in Mass 

Figure [3] compares cluster- forming cloud fragments to 
clouds devoid of significant clusters. For Taurus and the 
Pipe Nebula, the mass-size data is presented as done 
before (e.g.. Fig. [T|). We stress again that the Taurus 
MAMBO data for small spatial scales does not cover the 
entire cloud. At r < 0.2 pc, the data only character- 
ize the conditions in regions devoid of significant stellar 
groups. For Perseus and Ophiuchus, we choose a dif- 
ferent plotting scheme. Here, we only plot the data for 
fragments containing the most massive fragment found at 
small radii (NGC1333 in Perseus, L1688 in Ophiuchus). 

We find that, at given radius, cluster-forming Ophi- 
uchus cloud fragments (i.e., towards L1688) are signifi- 
cantly more massive than those in Taurus and Pipe. The 
naive expectation (i.e., that cluster-forming fragments 
are more massive) is thus confirmed. For Perseus (i.e., 
towards NGC1333), however, the result is more nuanced. 
At r < 0.3 pc, cluster-forming fragments are significantly 
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more massive than any structure in Taurus or the Pipe 
Nebula. For 0.4 < r/pc < 2, however, the extinction- 
derived mass towards NGC1333 is similar to the maxi- 
mum extinction-based masses for Taurus and the Pipe 
Nebula. If true, this would suggest that some cluster- 
forming regions at r > 0.4 pc have a structure similar to 
clouds only forming isolated stars (or no stars at all). 

It may be, though, that the Perseus extinction-based 
masses are biased towards lower values. The Perseus 
map has relatively poor physical resolution (5' at 260 pc 
distance), and intrinsic colors of stars in NGC1333 can 
bias extinction measurements. These problems are am- 
plified by high column densities towards NGC1333. Sec- 
tion 12.2.11 demonstrates these problems for B59. Thus 
we speculate that, at given size, the region containing 
NGC1333 is indeed more massive than any structure in 
Taurus and the Pipe Nebula. This remains to be proven, 
though. 

The fragments in the Orion Nebula Cluster, as well as 
Orion A as an entire cloud are more massive than any 
feature in Taurus and the Pipe Nebula (Fig. [2]). These 
data are thus consistent with the aforementioned trend. 

In summary, fragments containing clusters appear to 
be more massive than all structure in clouds devoid of 
clusters. Larger samples, and better data, are needed to 
ultimately establish this trend. 

3.4. Global Slopes of cluster-forming Regions 

As we just have seen in Sec. I3.3.1[ cluster-forming re- 
gions dominate the mass reservoir of their host cloud at 
any given spatial scale. This suggests a tight correlation 
between the properties of the cluster-forming fragments 
and the large-scale cloud structure. This can, for ex- 
ample, be characterized by "global" slopes of connection 
lines between mass-size measurements of the most mas- 
sive fragments at large and small scale (Eq. S]). Since 
we lack comprehensive data on all spatial scales, we are 
unfortunately presently unable to derive global slopes for 
GIO and Gil. To do so, one must also carefully charac- 
terize flux losses due to spatial filtering, as they occur in 
bolometer and interferometer maps. This is beyond the 
scope of the present paper. 

Specifically, the cluster-forming clouds are probed both 
at ~ 0.05 pc and ^ 5.0 pc radius (e.g., Fig.[T]). We calcu- 
late their global slopes by connecting the most massive 
fragments detected within 10% of these reference radii 
(30% radius deviation for the sparsely sampled Orion 
data). For Orion, where only very large scales are probed 
reliably, we use a mass measurement at 11 pc radius for 
slope calculations. Also, in Orion, the ^ 100 M0 frag- 
ments at 0.05 pc radius (right in the center of the Orion 
Nebula) are rejected as outliers. This rejection of objects 
is a regrettable move, since cloud fragments with extreme 
properties could be the actual sites where the most mas- 
sive stars are born. However, we feel that a more detailed 
and careful analysis than possible here, including a de- 
tailed consideration of the temperature structure, is war- 
ranted. We hope to do this in one of our future studies. 
The resulting slopes range from 1.10 to 1.33, as shown 
in Fig. 01 Note that a much smaller slope (as small as 
~ 0.7) would hold for Orion A, if we do not reject the 
outliers. The uncertainties indicated in Fig.Uhold for an 
uncertainty of a factor 2 in the mass ratio (factor 4 for 
the less reliable Orion data). This error budget presents 
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Fig. 4. — Global slopes for cluster- forming clouds. The slopes 
refer to the trends indicated in Fig.[T] the calculation of uncertain- 
ties is discussed in the text. The uncertainty- weighted mean slope 
of the sample is 1.27 ± 0.11. All clouds are consistent with having 
this slope, but the individual slope uncertainties are significant. If 
we include (highly uncert ain) Orion Nebula data, slopes as low as 
~ 0.7 are found. See Sec. l3.4lfor details. 



the worst expected scenario; the true errors are probably 
smaller. In Ophiuchus, we have sufficient resolution to 
also measure masses at ~ 0.02 pc radius. If we use such 
higher resolution data for this cloud, the slope decreases 
to 1.30. 

Note that we experimented with different approaches 
to define global slopes to characterize the core-cloud 
structure. We then decided to use the current method 
using masses only at 0.05 pc and 5.0 pc because we find 
no significant differences compared with the other meth- 
ods. For example, following the construction of Eq. 
we measure for each cloud the slope, b, that minimizes 
max[m(r)/(r/pc)'']. This value is then derived by uti- 
lizing data at all spatial scales. Slopes derived in this 
fashion only differ by ±0.04 from those shown in Fig. 
01 though. In the end, Eq. appears to provide the 
simplest approach that can also be repeated by other 
researchers. 

The uncertainty-weighted mean global slope for the 
three regions is 1.27 ± 0.11. As shown in Fig. [H our 
data are consistent with the hypothesis that all cluster- 
forming clouds have the same global slope. This sug- 
gests that mass-size relations m oc 7-i-27±o.ii pj-gvide a 
crude tool to use mass measurements at large spatial 
scale (^ 5 pc) to gauge the star formation conditions of 
the most massive embedded regions (~ 0.05 pc). Un- 
fortunately, the extreme SCUBA cores in Orion A, with 
masses ^ 100 Mq, suggest that much smaller slopes may 
hold in some clouds. The data are incompatible with 
larger slopes, though. Thus, the uncertainties shown in 
Fig. m suggest that the slopes are not larger than 1.5. 

3.5. Tangential Slopes: A Transition to Dense Cores? 

Figure |S1 presents tangential slopes for the entire sam- 
ple. Because of too strong spatial filtering, we exclude 
maps taken with interferometers, SCUBA, and Bolo- 
cam from this calculation. For MAMBO, Taurus frag- 
men ts larger r = 0.1 pc might be affected by filter- 
ing (jKauffmann et aLll2008D . Their slopes are thus not 
drawn. 

A first trend noted in Fig. |S1 is that the infinitesimal 
slopes change with radius. Also, they have no obvious 
relation to the global slopes shown in Fig. |4l Both slope 
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Fig. 5. — Tangential slopes for the sample clouds. See Fig.[T]for 
mark-up and coloring. The dotted line indicates the up per slope 
limit inherent to our measurement technique (Sec. [2rL2| l. Because 
of spatial filtering, slopes from Taurus MAMBO data may be bi- 
ased for r > 0.1 pc; such data are removed here. Data from other 
bolometers are not drawn because of the same problem. The black 
line gives the slope for B68, a Bonnor-Ebert-like dense core; similar 
slope trends are expe cted for other cores resembling Bonnor-Ebert 
spheres. See Sec. 13 Si for details. 

definitions are tlius complementary. 

At given spatial scale > 0.5 pc (in this domain, we 
only have data for Perseus, Taurus, and Ophiuchus), all 
clouds have comparable infinitesimal slopes. Those for 
Ophiuchus are a bit lower than for the other regions, in- 
dicating that the observed clouds differ in their structure 
at r > 0.5 pc. In all clouds the slopes decrease with in- 
creasing radius. As we detail in Sec. I4.3[ this behavior 
resembles the mass-size trends for model clouds with fi- 
nite extent (i.e., density vanishes outside a finite radius). 

In Taurus and the B59 and B68 regions of Ophiuchus, 
we derive very small infinitesimal slopes for r < 0.1 pc. 
(These regions are the only ones well probed at such 
small scales.) Near r = 0.1 pc, some fragments are 
observed to have slopes as low as d ln(m)/d ln(r) ~ 
1. This result is not entirely unexpected. Spheres 
supported by isothermal pressure, for example, obey 
d ln(TO)/d ln(r) = 1 at intermediate radii (Sec. I4.3p . 
Such (Bonnor-Ebcrt) spheres are commonly believed to 
constitute good idealizations of dense core structure. In 
particular, B68 is today quoted as the tex tbook model of 
a Bonnor-Ebert sphere (|Alves et al.ll200l"l ). As expected, 
B68 shows a slope approaching 1 for large radii. We cau- 
tion that the spatial filtering of bolometer maps, as well 
as the removal of extended features in extinction maps, 
might bias the slope measurements towards lower values. 
However, at least for cores with simple geometries, like 
B68, slopes ^ 1 are expected for larger radii. This is 
required for Bonnor-Ebert-like density profiles. If such 
slopes were not obser ved, the paradigm of Bon nor-Ebert- 
like dense cores (e.g., iBergin k, Tafallal 120071 ) would not 
be consistent with our data. 

Note that the Taurus and Perseus infinitesimal slopes 
at r ~ 0.5 pc are much larger than the slopes expected 
for the embedded dense cores of simple geometry (~ 1 
for every core at some radius). If this dichotomy is real, 
it has interesting implications. Basically, the slope must 



decrease from > 1.5 at r > 0.5 pc to d ln(m)/d ln(r) ~ 1 
for radii 0.1 pc. In other terms, there appears to be 
a transition in slope from the diffuse cloud structure to 
self-gravitating dense cores, when considering decreasing 
radii. Again we must caution that these trends might 
partially refiect artifacts in the data. As demonstrated 
before (Sec. 12. 2H) . the mass in extinction maps is some- 
times increasingly underestimated towards small scales. 
This biases slopes towards lower values. It is not likely, 
though, that all features in our extinction maps are mas- 
sively biased. In particular for Perseus, such would be 
inconsisten t with comparisons be tween CO and extinc- 
tion maps (jGoodman et al.l [20091 : also see part I). Still, 
it is discomforting that the slope transition does appar- 
ently exactly occur in the spatial domain that is not well 
probed by any tracer. 

Interestingly, a cloud-to-core slope transition might not 
exist in Ophiuchus: note that slopes ~ 1 are observed for 
radii > 0.3pc, so that no slope transition between cloud 
and core might be necessary. Specifically, at any given ra- 
dius, the most massive fragment containing L1688 follows 
tightly the line connecting the most massive fragments 
at 0.05 pc and 5 pc radius. If this is true, a clean divi- 
sion between the L1688 dense cores and the surrounding 
diffuse cloud structure might not exist. 

4. INTERPRETATION OF MASS-SIZE DATA 

4.1. Larson's 3^'^ Law and other previous Work 

iLarsonl (jl981h carried out one of the first studies of 
mass-size relations of molecular clouds (in his Fig. 5) and 
obtained m{r) oc r^'^. (Actually, he studied size-density 
relations. Since he assumed spherical cloud geometries, 
these can be turned into mass-size laws.) Today, this 
"third Larson law" is usually quoted as m(r) oc r^ and 
considered one of the funda mental properties of cloud 
structure (e.g.. lMcKee fc Ost rikcr 200l). 

We find a nuanced relation between our data and the 
detailed result of Larson (.19811 ). 

m(r) = 460 M0(r/pc)i-^ 

(our Eq. [T]). In the 1 < r/pc < 4 radius range, his and 
our results are compatible for several solar neighborhood 
clouds, as shown by Eq. ([7]). However, some clouds, like 
GIO, can deviate by about an order of magnitude in mass 
from this relation. Further, some trends within individ- 
ual clouds are not described well by one single mass-size 
law. Slope and intercept of Eq. ([6|), and the slopes shown 
in Figs. [3] and [5l m ay serve as exa mples for deviant mass- 
size relations. The Larson! ([1981') mass-size law thus fails 
to describe a lot of the cloud substructure visible in state- 
of-the-art cloud maps. Devia tions from this relation h ave 
also been noted before (e.g., |McKee fc Ostrikeill2007l for 
references). Also, note that Larson (1981) never intended 
to match all cloud substructure; as seen in his Fig. 5, the 
mass-size relation only provides an order-of-magnitude 
fit to his data. In this sense, the mass-size laws derived 
here do not really supersede Larson's results; molecular 
clouds are very complex structures, and every one has to 
be considered individually. However, for certain mass- 
size trends, Eq. ([6|) and the global slopes in Fig. |4] do 
provide a more detailed description than offered by Eq. 

Beyond these quantitative differences, Larson's tech- 
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nique differs fundamentally from ours. He plotted all 
cloud and dense core data available to him, and at- 
tempted to fit this ensemble data by a common law, in- 
dependent of whether the cores and clouds were residing 
in the same area of the sky. Here, however, we usually 
handle data separately for each cloud, to e.g. calculate 
the global and tangential slopes shown in Figs. |3]and[5l 
Only Eqs. ([71 16]) stem from an analysis of data combined 
for several clouds. 

Ensemble-based mass-size relations of the kind con- 
sidered by Larson have been studied by many other 
groups ; here we ca.n on ly present excerpts from such 
work. iLada et"aLl ()2008[ ). for example, derive b = 2.6 
for dense cores in the Pipe Nebula extinction map used 
here. For a sampl e of Infrared Dark Clouds (IRDCs), 
iRagan et al.l (|2009f ) derive even steeper slopes, 6 w 3, 
much in excess of any slope derived here. Their ap- 
proach is very different from ours, though; they extract 
dense cores, plot their masses and sizes, and fit the data 
with a single line. We, instead, measure a core's mass 
at various contours, and derive tangential slopes by only 
considering data for a single dense core. Their ensem- 
ble slopes thus serve a different purpose than our global 
and tangential slopes. In this respect, it may be help- 
ful to consider the Ty pe 1-4 linewidth-size relations of 
iGoodman et al.l ()1998l ; replace line width with mass in 
their discussion): there are many different ways to de- 
fine slopes, and they will characterize different proper- 
ties. The large-scale structur e of m olecular clouds, as 
originally addressed bvlLarsoi] (|1981l). is in princi ple dis- 
cussed bv lSolomon et al.l (|1987D and lHever et al.l (|2009[ ). 
Their data has, however, not been exploited to derive 
mass-size rel ations. 

Note that iStuwd (|1990l ) uses a technique similar to 
ours to probe the spatial domain considered here. As 
discussed in Sec. 14.31 his work is consistent with ours. 

4.2. Mass-Size Relations and Star Formation 

It is thought since long that clouds need to achieve a 
high c olumn density in order to produce dense cores and 
stars. iJohnstone et al.l (|2004D . in particular, introduced 
the concept of extinction (or column den sity) thresh- 
olds for dense core formation (also see: [Onishi et al.l 
iggs", 'HatchcU ct al." 2005, Enoch ct al. 2007). Similarly, 
Lpmbardi ct al. (2006) and Lada et al. (2009) argue that 
a low fraction of mass at high colu mn density yields 
low s tar formation activity (also see iKainulainen et all 
l2009t) . This is in line with the Kennicutt- Schmidt law 
between star formation rate and mass surface density, 
SsFR oc (e-g-, Kennicutt 1998). Since p ^ 1 for star- 
forming clouds ((Evans et al.il2009[ ). this relation predicts 
an increase of star formation activity with increasing col- 
umn density. 

The analysis in Sec. [3] shows that these laws do also 
manifest in our data. Sections l3.ll and [331 suggest that 
the ability to form clusters and massive stars increases 
with increasing mass, when considering a given radius. 
Since {NH2) oc m/r^, one could also say that cloud frag- 
ments do appear to only form massive stars and clusters 
when they have a high mean column density — ^just as 
suggested by the aforementioned laws. 

Mass-size studies do, however, also provide informa- 
tion not available from the aforementioned plain column 



density studi es. First, note that "e xtinction threshold" 
studies (e.g., IJohnstone et al.l [20041) do only consider a 
single spatial scale, i.e. the beam used to construct the 
extinction map. This is a major difference to mass-size 
studies, where many spatial scales are considered. Sec- 
ond, observe that studies of colum n density distributions 
(PDFs; e.g.. lLombardi et al]|2006[) do consider all scales 
of a map, but do not register which part of the signal 
shown on the histogram originates is which part of the 
cloud. (Depending on the analysis, this is not necessarily 
a problem.) Mass-size studies, in contrast, treat individ- 
ual cloud fragments separately. 

Mass-size data sets give a new twist to discussions 
of extinction thresholds. To see this, consider the so- 
lar neighborhood clouds examined here. These are all 
similar at large spatial scale (Eq. [7|). Still, at smaller 
scale, they differ significantly in mass and star formation 
activity (Sec. 13. 3|) . The processes determining the star 
formation activity must thus operate on spatial scales 
smaller than the entire cloud. In this sense, the star for- 
mation activity depends (at least in our sample clouds) on 
a cloud's ability to create, from a given mass reservoir, a 
small number of fragments that dominate the mass reser- 
voir and concentrate it into a small volume. The pres- 
ence of high column densities are then a consequence of 
the cloud structure, but not the governing reason for the 
formation of dense cores and stars. 

4.3. Slopes and Intercepts: Constraints on Density and 
Physical Cloud Models 

The observed power-law-like mass-size relations, 

m(r) = tuq (r/pc)^ , 

are characterized by slopes, &, and intercepts, tuq. As we 
show here, slopes and intercepts can be used to gauge 
densities and their gradients. At the same time, it is 
possible to constrain the absolute level of pressure, and 
the nature of its origin. 

Consider, for example, an infinite equilibrium sphere 
with power-law density profile that is supported by 
isothermal pressure from gas at temperature Tg. Then, 

™M = 2.6Mo(i)(5^) (8) 

(jKauffmann et al.ll2008l Eq. [13], in their case e 7r/2). 
Thus, if this model holds, the intercept encodes informa- 
tion on the gas temperature supporting the cloud. Con- 
versely, the slope predicted by the model, 6=1, can be 
used to validate the model; if observations yield 6 7^ 1, 
then the model does not apply. 

It has to be kept in mind that we consider mass-size 
laws derived from two-dimensional maps. These are 
related to, but not identical with, mass-size laws ob- 
tained from three-dimensional density map s. This is il- 
l ustra ted by the experiments conducted bv lShettv et al.l 
|2010() who use the fragment identification technique also 
used by us. Their analysis is based on three-dimensional 
numerical simulations of turbulent clouds. As part of 
their experiments, they derive mass-size slopes from their 
data. For their particular set of simulations, the power- 
law slopes derived in this fashion are similar to the num- 
ber of dimensions used for mass measurements (i.e., 3 
when based on density, and 2 when based on column 
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density). This underlines that the number of dimensions 
considered has to be kept in mind. 

4.3.1. Density Laws 

The above discussion can be extended to include many 
more models of cloud structure. We do this in two Ap- 
pendices. The results of this analysis are summarized 
here. 

First, let us examine the connection between mass-size 
laws and cloud density profiles. Consider a sphere with 
about constant density for radii smaller than some flat- 
tening radius, sq, a power-law drop at intermediate radii, 
n{s) oc s~'' (where s is the distance from the center), and 
vanishing density beyond some outer truncation radius, 
R. Such profil es provide a good match to observations 
of dense cores (|Tafalla et al.il200l) . They are a good ap- 
proxima tion to the structur e of isothermal equilibrium 
spheres (jPapp &: Basull2009( ). As we show in Appendix 
lA.H for apertures of radius r, this yields mass-size rela- 
tions of the form 

m{r) oc Uc r^"*" for sq < J' < ^ , (9) 

where ric is the density for s = 0, and mass-size slopes 

d In(TO) /d ln(r) = 3 - k for sq ^ r <^ R . (10) 

Both relations apply only if fc < 3. As the equations 
show, the intercept contains information on the central 
density, and the mass-size slope depends on the slope of 
the density law, k. Both does, of course, only hold at 
intermediate radii, sq <^ r <^ R. For reference, we note 
that the column density obeys N oc r^~''. 

A generalized versio n of p ower-law spheres are tri-axial 
ellipsoids. Appendix IA.3I considers the case in which 
n(s) oc (s/sq) ^^ along any main a xis, but with sq de- 
pending on the direction chosen (Eq. IA9I) . Detailed anal- 
ysis shows that such ellipsoids follow the same mass-size 
relations as spheres, when r = (A/ny/^. Thus, the laws 
listed above apply. 

In a next step, one may wish to consider models of 
cylindrical clouds of length £. Here, we adopt density 
drops n{s) oc perpendicular to the cylinder axis for 
intermediate values of s. At intermediate radii so ^ r <C 
R, such clouds obey 



m(r) oc ric r 



,4-2fc 



d ln(m)/d ln(r) = 4 - 2fc (11) 



if their major axis is perpendicular to the line of sight, 
and 



m{r) oc Tic r 



d ln(m)/d ln(r) 



(12) 



if the axes are aligned. (Intermediate angles are not con- 
sidered here.) Meaningful relations are only obtained 
for k < 2. In both relations, the radius is defined as 
r ^ (A/7r)i/2. 

Thus, just as one may naively expect, the mass-size 
slope gauges the slope of the density profile. Further, 
the intercepts of mass-size relations constrain the abso- 
lute density of cloud fragments. There is, however, one 
less obvious fact that calls for attention: the exact re- 
lations between mass-size slopes, intercepts, and density 
law depends on the cloud model and viewing angle. It is 
therefore not possible to derive the true density profile 
without further information on the cloud geometry. Such 



information may, e.g., be derived by studying the elon- 
gation of cloud fragments. Also, the above power-law re- 
lations do only apply at intermediate radii, sq <C r <C i?. 
This domain might not exist in actual observed clouds. 
Then, the central density plateau and the finite size have 
to be considered. These give 



m(r) 



oc ric J' 
w M 



^ for r ^ So and 
for r > R . 



(13) 



4.3.2. Polytropic Equilibria 

The density slopes themselves depend on the processes 
shaping the model cloud. As a first example, here we con- 
sider static equilibrium models in which pressure gradi- 
ents are in balance with self-gravity. We assume a poly- 
trophic equation of state, P oc n^^, in which pressure 
and density are related by the poly trophic exponent, 7p. 

In Appendix |B] we show that 



k = 



2-7P 



(14) 



for polytropic equilibrium spheres (if jp < 4/3) and 
cylinders (if 7p < 1). The density and mass-size slopes 
are, thus, related to the polytrophic exponent. Isother- 
mal pressure, for which 7^. = 1, implies fc = 2 in spheres, 
for example. Then, d ln(m)/d ln(r) = 1 in spherical 
model clouds; laws too complex to be considered here 
apply to cylinders. As seen in Figs. |4] and [SJ such a model 
can explain some, but not most slope measurements. 

Polytropic exponents 7p — 1/2 are sometimes sug- 
gested to describe "turbulent" pressure within clou ds, as 
e.g. arising from Alfven waves (IMcKee fc Zweibeli riOgSl. 
In this case, fc = 4/3, and so d ln(?7i)/d ln(r) assumes 
values of 5/3 ~ 1.67 (spheres and ellipsoids), 4/3 « 1.33 
(perpendicularly viewed cylinder), and 2/3 ~ 0.67 (end- 
on cylinder) for the different models. Among these, 
spheres, ellipsoids, and cylinders viewed from the side 
provide an acceptable match to the observed mass-size 
slopes b > 1. Cylinders viewed along their major axis 
yield too shallow mass-size laws (and such a viewing di- 
rection is highly unlikely). 

For a given physical model, the intercept can be used 
to gauge a cloud's stability against collapse, respectively 
suggest the level of supporting pressure. Here, we limit 
ourselves to the isothermal case, i.e. = 1. Stability 
consideration s (e.g., of Bonnor-Ebert-type; lEbertl [19551 
lBonnoHri956D imply 



M < Mr, 



a^{v)R 
G 



(15) 



for the total mass, where a{v) is the characteristic one- 
dimensional velocity dispersion fEqs. [B3l and iBS]) . For 
spheres, R is the radius, while R ^ £ in cylinders. If a{v) 
is known (e.g., for thermal pressure), M > A/cr implies 
collapse of the object considered. Conversely, depending 
on the situation, a{v) can be inferred by requiring that 
AI = Met- Required values of a{v) significantly in ex- 
cess of the thermal velocity dispersion of the mean free 
particle might, e.g., suggest the presence of significant 
non-thermal pressure. If we only require that pressure 
balances gravity, and drop the constraint that the object 
is to be stable against perturbations, the above law yields 
Eq. (HI). 
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As particular example, consider B68. It is thought that 
this dense cor e has a structure very similar to a Bonnor- 
Ebert sphere (|Alves et al.ll2001[ ). Thus, one would expect 
the mass and size of B68 to obey Eq. ([8|) , when consid- 
ering large enough radii. This is indeed the case, as seen 
in Fig. m 

4.4. Synoptic and physical Density Slopes 

For intuitive communication, it may be helpful to re- 
port synoptic density slopes. 



d ln(n) 



d ln(s) 



d In(TO) 
d ln(r) 



(16) 



i.e. the density slope a sphere of the same mass-size slope 
would have when observed at intermediate radii. The 
synoptic slopes give a good first idea of the true density 
slopes. First, recall that the model mass-size laws do not 
sensitively depend on the assumption of exact spheres; 
the same relation holds for ellipsoids. Also, in the ob- 
served range 1 ^ d ln(TO)/d ln(r) < 2, spheres (or ellip- 
soids) and perpendicularly viewed cylinders (the end-on 
view is statistically insignificant) imply similar slopes; 
for these geometries, the synoptic slopes exceed the true 
ones by less than 0.5, assuming intermediate radii. Thus, 
we derive 



-d ln(n)/d ln(s)]^. 



syn 



1 to 2 



for mass-size slopes 1 < d ln(m)/d ln(r) < 2. 

A limited comparison of these de nsity-siz e slope s with 
previous results is possible. Tafall a~et al.l (|2002l ). e.g., 
study the dust emission of five starless dense cores, and 
derive density-size slopes of 2.0 to 2.5 for four of them. 
This is a typical result for c loud fragments of < 0.1 pc 
size (iBergin fc Tafallal 120071 ). a spatial domain not too 
well covered by our data. 

Concerning the analysis meth od and spatia l range con- 
sidered, the extinction study bv lStiiwd ()1990[ ) might pro- 
vide the best match to our work. Based on star counts, 
he derives — d ln(rt)/d ln(s) > 1.0 ± 0.4 on scales of up 
to ~ 1 pc. This is consistent with our results, also given 
the differences in map construction (he uses optical star 
counts). 

5. SUMMARY & OUTLOOK 

This work studies the internal structure of molecular 
clouds by breaking individual cloud complexes up into 
several nested fragments. For these, we derive masses 
and sizes in order to study their density structure. 

Analysis of a limited sample of solar neighborhood 
cloud complexes (Taurus, Ophiuchus, Perseus, Pipe Neb- 
ula, Orion A), as well as more distant clouds (GIO & 
Gil), yields first some basic constraints on mass-size 
cloud structure. These are as follows. 

1. On large spatial scales, > 1 pc, the most massive 
fragments in solar neighborhood clouds — with the 
possible exception of Orion A — obey 

mir) = 400 Mq (r/pc)^-^ 

(Eq. [7]) with deviations < 40%. This relation 
resemb le s the original mass-size law derived by 
ILarsonl (pSl ). m(r) = 460 Mq (r/pc)i ^ The 



more distant clouds in the sample, however, deviate 
from this relation by up to an order of magnitude 
in mass. 

2. No single mass-size law can be used to describe all 
fragments in all clouds. In particular, "Larson's 3'''^ 
Law" of constant column density, m(r) oc r^, pro- 
vides a bad global description; today's data are too 
complex to warrant the use of such relations. To 
give examples, power-law slopes vary with radius 
within a given cloud (Fig. [5]), and clouds can differ 
massively in mass at given radius (Fig. [2]). 

In practice, different definitions of mass-size laws 
are used by different researchers. Also, different 
definitions may serve different purposes. This must 
be taken into account when comparing different 
mass-size laws. 

Most importantly, the mass-size data can be used to learn 
about the formation of stars in molecular clouds. We 
derive the following constraints. 

3. Sample clouds not forming massive stars (> 
10 Mq) adhere to a limiting mass size relation. 



m(r) < 870 Mq {r/^c 



1.33 



(Eq.|ni), while our sample clouds forming such stars 
violate this law (Fig. [5]). This suggests that the 
above relation describes the typical mass-size range 
of molecular clouds not forming high-mass stars. 
Also, the observations advocate that this bound- 
ary constitutes a mass limit for massive star for- 
mation. However, such conclusions are based on a 
small sample and are thus preliminary. 

4. Across all clouds studied here, cloud fragments 
forming clusters are more massive than fragments 
not doing so (Figs. [T] and El Sec. 13.3.1]) . At given 
size, cluster-forming fragments dominate the mass 
reservoir of their host cloud. 

5. The mass-size trend of cluster- forming fragments 
can e.g. be captured by global mass-size slopes 
(i.e., from 0.05 pc to 5.0 pc radius; Sec. 13. 4[) . Our 
cluster-forming sample clouds are consistent with 
a common slope ^ 1.27. The uncertainties are, 
unfortunately, significant for a given cloud; slopes 
may well differ between clouds. In the case of 
Orion A, e.g., the slope might be as low as ^ 0.7. 
However, in any event slopes smaller 1.5 do hold. 

Theoretical discussions show that mass-size laws of the 
form m{r) = mo (r/pc)^ can be related to physical cloud 
models characterized by power-law density gradients, 
n{s) oc s~^, or polytropic equations of state, P oc n^'' 
(Sec. 14.31) . Provided certain idealizations apply, 6, fc, 
and 7p are directly related to another. This analy- 
sis suggests the definition of a synoptic density slope, 
[— d ln(n)/d ln(s)]gy„ = 3 — d ln(m)/d ln(r) (i.e., assum- 
ing the fragment considered was a sphere). This slope 
provides a first rough estimate of the true density law. 
Our data gives synoptic density slopes in the range 1 to 2. 

We are indebted to a careful and thorough anonymous 
referee, who saved us from making an embarrassing mis- 
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take. This project would not have been possible with- 
out h elp from Erik Rosolowsky . His dendrogram analysis 
code (iRosolo wskv et al] l2008t) was instrumental for our 
analysis . Carey et al. (I2000D. Alves e t al. (priv- comm.) 
lEnoch et a l. (2006), 'Lombardi ct al.' (2006VR.idge_eLalJ 
(|2006( ). iRoman-Zuniga et al.. (2009. ). and Thompson et 



al. (in prep.) contributed maps to the present study. We 
are grateful for their help. This work was in part made 
possible through Harvard Interfaculty Initiative funding 
to the Harvard Initiative in Innovative Computing (HC). 
It is based upon work supported by the National Science 
Foundation under Grant No. AST-0908159. 
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APPENDIX 

A. MODEL MASS-SIZE RELATIONS 
Here we consider density profiles of the form 

n{s) oc s^'' (Al) 



in the domain sq ^ s R. In spheres, s is the distance from the density peak, while it is the distance from the main 
axis in cylinders. Densities become about constant for s < sq, and they vanish for s > R. 
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A.l. Spheres at Intermediate Radii 
Integration along the line of sight gives the column density at offset u, 

N{u) = 2 / n{[u^ + w^Y'^)dw , 

Jo 

where the integration stops at the boundary radius. Evaluation yields 



(A2) 



(A3) 



if A; > 1, a relation that can be used for offsets sq ^ u R. Integration of the column density over a circular 
peak-centered aperture of maximum offset u gives the mass, 



m{u) ^ 2ttC / N{w)wdw. 
Jo 



(A4) 



where C = Yi/Nn^ — /ih2"^h is the conversion factor to mass surface density. Evaluation in the limit i? — > oo gives 

m(u) cx u^^'^ m(r) cx r^^^ (A5) 

for sq <^ r R. We use u = r in the transition from offsets to aperture sizes since, in spheres, the offset from the 
center, u, is equal to the radius defined via the aperture area, r = {A/ttY^^. 

The meaningful range of this relation is limited to fc < 3; for fc > 3, the mass would decrease with size, which is not 
physical. In essence, spheres with A: > 3 have a finite mass, even for R — >■ c3o, so that power-law mass-size relations 
cannot hold in any radius domain. 

A. 2. Homogeneous Spheres 

Homogeneous spheres, in which the density is spatially constant, constitute a special case of the aforementioned 
spherical models. Consider an offset u from the sphere's center, where the density drops to zero beyond the sphere's 
outer radius, R. Then, the column density is the product of the density, no, and the length of the line of sight. 



Integration following Eq. 



4]) yields (with the substitution u 
4 



m{r) 



r3 _ (i?2 _ ^2)3/2 



(A6) 



(A7) 



in which go = Cuo is the mass density corresponding to no- For r ^ R we derive m = 4/37r£io^^, just as expected for 
a homogeneous sphere truncated at radius R. The slope obeys 



d ln(m) 
d ln(r) 



(i?2 _ r^) 



2U/2 ^2 



i?3 _ (i?2 _ ^2)3/2 



(A8) 



and does monotonously decrease from d ln(m) /d ln(r) = 2 at r = to d ln(m) /d ln(r) — at r — R. 



A. 3. Triaxial Ellipsoids at Intermediate Radii 
Another variation of spherical power-law density profiles are triaxial density distributions with ellipsoidal iso-density 
surfaces. These can be derived from spherical power-laws by substituting s/sq — [(x/xq)^ + {y/yoY + '^oYY^'^ 
(where a subscript '0' indicates reference properties for normalization). The coordinates x, y, and z give the projection 
of a given position on the three axes of the ellipsoidal density distribution. Then, the density law reads 



no 



-k/2 



(A9) 



The chosen coordinates form an orthogonal coordinate system. Any position along a straight line of sight, w, is 
therefore linearly related to the coordinates chosen for the density distribution: 



fTT-x + ; y = niy + UyW; z = niz + UzW . 



(AlO) 



For c onven ience, her e we choose that w = at the position of highest density along the line of sight. Substitution of 
Eqs. (KTO\i into Eq. yields 



n = Ho 



ml ml ml 
^0 Va ^0 



m^Ux myUy m^n^ 



Xo 



yo 



zo 



-k/2 



(All) 
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i.e., some terms that do not contain w, some that include w hnearly, and a few that contain w to its square. Analysis 
shows that the sum of the terms with a linear dependence on w must be zero. To see this, consider the requirement 
that the density reaches its maximum where w = 0. For this to happen, the derivative of the sum within the square 
brackets o f Eq. ETD with respect to w must vanish at w = 0. This is only the case if all terms linear in w cancel out. 
Equation (jAlip can thus be written 

n = no ( — j , (A12) 

if one chooses = (m^/xg + ml/y^ + ml/z^) ■ {nl/xl + n^/yg + nl/ and rg = {nl/xl + nl/yl + nl/zly^. 
Formally, Eq. (jA12l) is identical to the spherical density law substituted in Eq. (jA2|) . Therefore, we can use the 
calulations in Appendix I A. 1 1 to conclude that, again, 

N{u) oc u^-^ {R oo) (A13) 

\i k > \. Thus, the spatial column density distribution does only depend on u. 

A detailed look at the above analysis does actually show that we did not use the exact nature of the density law 
to this point. Instead, we exploited that the density depends on (x/xq)^ + (y/j/o)^ + (-z/^^o)^- Thus, we can make the 
conclusion that all density distributions with ellipsoidal iso-density surfaces yield column density distributions only 
depending on u. 

Thus, u can be interpreted as an offset, just as in the spherical situation. Inspection of the equation defining u 
reveals that the three-dimensional offsets from the distribution's center, m^, form a triaxial ellipsoid too, when u is 
kept constant. Projection of the u — const, surface into a plane corresponds to a cut through this ellipsoid. Since cuts 
though ellipsoids yield ellipses, the u = const, curve in a plane is an ellipse too. In this case, spatial integration of the 
column density can be conveniently executed in polar coordinates. Let v and w be the minor and major axis of the 
projected ellipse. For a given viewing perspective, v /w is constant for all contours of constant column density. Then, 
the mass follows from ^ 

m{u) ^2ti (v/w)C j N{w)wdw, (A14) 
Jo 

where we use that the ellipse's area increases with w as 27r (v/w) w dw. This integral becomes 

m(u) oc u^^'' =^ m{r) oc r^~^ (A15) 

for sq ^ r R. In this, we use that u is proportional to the effective radius of the ellipsoid formed by iV = const, (i.e., 
u — const.) contours. The latter holds since, for a given ellipsoid, v/w is a constant, and its area is -kvw = tt{v /w)w'^ . 
Thus, the effective radius becomes r = {v/wY^^w oc w. 

In summary, we retrieve the relations already derived for spheres, including the limits on reasonable values oi k. As 
expected, spheres form a particular case of ellipsoids with xq — uq — zq. 

A. 4. Cylinders at Intermediate Radii: Perpendicular View 

In cylinders with the main axis perpendicular to the line of sight, the integration in Eq. (|A2p must be evaluated 
with the integration path perpendicular to the cylinder axis. Then 

N{u) oc u^-'' (R^oo), (A16) 

as long as fc > 1. This provides, again, a convenient description for sq <^ u <^ R. Since the column density is constant 
for given offset, u, apertures along a given column density contour have a size 2u£, where £ is the cylinder length. Such 
an aperture contains a mass 



m{u) = 2£C N{w)dw. (A17) 
Jo 

Substitution yields 

m{u) oc u^^'' ^ m(r) oc r'^^^'' . (A18) 

This relation is applicable if fc < 2, following the discussion of finite masses for spheres. In the transition from u to r 
we use that r = {2£u/TTy^^, which gives u (x r^. 

A. 5. Cylinders at Intermediate Radii: Parallel View 

In cylinders with the main axis aligned with the line of sight, the integration in Eq. (|A2[) must be evaluated along 
the cylinder axis. In this case, 

N{u)^n{u)-£. (A19) 

This geometry further implies 

m{u) = 2TiC [ N{w)wdw, (A20) 
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which gives 

m{u) cx u'^^'^ => m{r) cx r^"'' (A21) 

in the usual sq <C u ^ i? limit, with the further constrain that k < 2. As in spheres, u = r for contours of constant 
column density. 

A. 6. Very large and small Radii 

For central flattening of the density profile, the column density becomes constant for u <C sqj independent of the 
model adopted. Then, the mass in a given aperture is just given by the product of the mass surface density, which 
scales with rici?, and the aperture area. For large radii, the mass is equal to the total mass when considering apertures 
larger than the object's size. Thus 

m(T\l '^^^^^^ for r<so , and 

"'^'^^ \ « M for r>R , ^^^^> 
independent of the adopted model geometry. 

B. POLYTROPIC EQUILIBRIUM CLOUD MODELS 
Polytropic cloud models assume an equation of state of the form 

P = P^ {Q/Qay" ■ (Bl) 

They are in particular used to solve the equation of hydrostatic equilibrium, VP = — p V^, where the gravitational 
potential fulfills the Poisson equation, V^<I> = AirGg. 

B.l. Polytropic Equilibrium Spheres 
Hydrostatic equilibrium spheres supported by polytropic pressure have recently been summarized by 
iMcKee fc HollimanI (|1999( ) and iCurrv fc McKed (|2000D . In spherical symmetry, V d/ds, $(s) = -Gm{s)/s, 
and m(s) = in g{w)wdw. This differential equation can easily be solved for density profiles g{s) oc s^'". Solutions 
are, however, only physical for fc < 3; the power-law implies m(s) oc s'^~^, which only increases with radius if k is 
small enough. The hydrostatic equation implies qc s^"^'^. This can only be fulfilled if both exponents to s 

are identical, which implies 

k = 7r^ (7P<4/3), (B2) 
2 - 7p 

where the limit on jp reflects the condition fc < 3. Stability ag ainst perturbations implies that the mass of the 
equilibrium is smaller than some critical mass (jMcKee fc HollimanI 1 1999,) . For 7^ = 1, this critical mass is identical 
to the Bonnor-Ebert mass, 

Mc. = 2.4^, (B3) 
where the numerical constant follows the discussion bv IMcKee fc HollimanI (|1999f ). Stable clouds have M < M„. 

B.2. Polytropic Equilibrium Cylinders 

iHoredtj (|1987l ) presents a summary of hydrostatic equilibrium cylii iders (as well as sheets and spheres) supported by 
polytropic pressure. It bui l ds on discussions of the case 7_p > 1 bv lOstrikeil (|1964[ ). and analysis of polytropes with 
7P < 1 bv lViala fc HoTedtl ([1971 . 

The Poisson equation can be tackled using the divergence theorem, JyVFdV = J^yFdA. When analyzing 
this equation for the volume of an infinitely long cylinder, the ends can be neglected m the integration over the 
volume's surface, dV. Thus, if the component of F perpendicular to the cylinder surface is constant, and assumes the 
value F± along it, then |/g^FdA| — 2tt£sF± for cylinders of length £ and radius s. To apply this analysis to the 

Poisson equation, we substitute F — >■ V<I>. In cylindric coordinates, the component of the potential's gradient that is 
perpendicular to the cylinder surface implies the substitution F± — d^/ds. In essence, combination of the Poisson 
equation and the divergence theorem yields AttG Jy gdV — Jgy d^/dsdA, which evaluates to d^/ds = 2Gm{s) / (si) , 

where m{s) is the mass within the radius s. For £ — 00, we further have V — > d/ds, and so the hydrostatic equation 
becomes dP/ds = —gd^/ds. 

A simple solution is provided by g{s) oc s~^ , which implies s~''^''^^^ oc s^^"^^ , just as found for the spherical case. 
Analysis of the integral for the mass implies fc < 2, similar to what is found for spheres. Thus, solutions must fulfill 

fc=-^ (7P<1), (B4) 
2 - 7p 
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where the hmit enforces fc < 2. A perturbation analysis analog to the spherical Bonnor-Ebert case gives a critical 
limiting mass 

Me,. = 2 — ^ (B5) 
for 7p = 1, as e.g. demonstrated by lOstrikeil (|1964[ ). 



